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above 2 GeV. That indicates that the Goldstone boson exchange effects may play 
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of [it<i] 2 c pentaquark ~ 3250 MeV and [ud] 2 6 pentaquark ~ 6509 MeV. 
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1. Introduction 

Recently LEPS and DIANA collaborations [U Ej reported the observa- 
tion of a very narrow peak in the K + n and K°p invariant mass distribution 
which existence has been confirmed by several experimental groups in var- 
ious reaction channels [3j. These experimental results were motivated by 
the pioneering paper on chiral soliton model |4J. The reported mass de- 
terminations for the are very consistent, falling in the range 1540±10, 
with the width smaller than the experimental resolution of 20 MeV for 
the photon and neutrino induced reactions and of 9 MeV for the ITEP 
K + Xe — > K°pXe' experiment. 

From the soliton point of view the is nothing exotic compared with 
other baryons - it is just a member of 10^ multiplet with S = +1. 
However, in the sense of the quark model + (154O) baryon with posi- 
tive amount of strangeness is manifestly exotic - its minimal configuration 
can not be satisfied by three quarks. The positive strangeness requires 
an s and qqqq (where q refers to the lightest quarks (u,d)) are required 
for the net baryon number, thus making a pentaquark uudds state as 
the minimal "valence" configuration. Later NA49 collaboration at CERN 
SPS announced evidence for an additional narrow uduss resonance 
with / = 3/2, a mass 1.862 ± 0.002 GeV and a width below the de- 
tector resolution of about 18 MeV 1 ) and HI collaboration at HERA P 
found a narrow resonance in D*~p and D* + p invariant mass combinations 
at 3.099 ± 0.033 sta t ± 0.005 syst GeV and a measured Gaussian width of 
12 ± 

3 s tat MeV, compatible with the experimental resolution. The later 
resonance is interpreted as an anti-charmed baryon with a minimal con- 
stituent quark composition of uuddc, together with the charge conjugate. 
The discoveries of first manifestly exotic hadrons mark the beginning of a 
new and rich spectroscopy in QCD and provide an opportunity to refine 

1 ^NA49 also reports evidence for a S°(1860) decaying into S(1320)7r. 
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our quantitative understanding of nonperturbative QCD at low energy. 

The O-hyperon has hypercharge Y = 2 and third component of isospin 
h = 0. The apparent absence of the I3 = +1, 6 ++ in K + p argues against 
7=1, therefore it is usually assumed the O to be an isosinglet. The other 
quantum numbers are not established yet. 

As to the theoretical predictions we are faced with a somewhat am- 
biguous situation, in which exotic baryons may have been discovered, but 
there are important controversies with theoretical predictions for masses 
of pentaquark states. The experimental results triggered a vigorous theo- 
retical activity and put a renewed urge in the need to understand of how 
baryon properties are obtained from QCD. 

All attempts of the theoretical estimations of the pentaquark masses can 
be subdivided into following four categories: (i) dynamical calculations us- 
ing the sum rules or lattice QCD [3, Ej , (h) the phenomenological analyses 
of the hyperfme splitting in quark model jHl HDj, (hi) phenomenological 
analyses of the SU (3) p mass relations, and (iv) dynamical calculations 
using the chiral SU(3) quark model [fTTj ) . 

The QCD sum rules predict a negative parity 6 + of mass ~ 1.5 GeV, 
while no positive parity state was found [7j. The lattice QCD study also 
predicts that the parity of the lowest G hyperon is most likely negative 

0. 

The naive quark models, in which all constituents are in a relative S- 
wave, naturally predict the ground state energy of a J p = | pentaquark 
to be lower than that of a J p = | + one. However, using the arguments 
based on both the Goldstone boson exchange between constituent quarks 
and color-magnetic exchange it was mentioned that the increase of hy- 
perfme energy in going from negative to positive parity states can be 
quite enough to compensate the orbital excitation energy ~ 200 MeV. 
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However, existing dynamical calculations of pentamasses using the chiral 
SU(3) quark model (see, e.g., [TT] ) are subject to significant uncertainties 
and can not be considered as conclusive. 

Pentaquark baryons are unexpectedly light. Indeed, a naive quark 
model with quark mass ~ 350 MeV predicts 6 + at about 350 x 5 = 1750 
MeV plus ~ 150 MeV for strangeness plus ~ 200 MeV for the P-wave 
excitation. A natural remedy would be to decrease the number of con- 
stituents. This leads one to consider dynamical clustering into subsystems 
of diquarks like [ud] 2 s [12] and/or triquarks like [udjfwds] (ID] which am- 
plify the attractive color-magnetic forces. In particular, in (12] it has been 
proposed that the systematics of exotic baryons can be explained by di- 
quark correlations. 

The quark constituent model have not been yet derived from QCD. 
Therefore it is tempting to consider the Effective Hamiltonian (EH) ap- 
proach in QCD (see, e.g., (13]) which from one side can be derived from 
QCD and from another side leads to the results for the qq mesons and 
3q baryons which are equivalent to the quark model ones with some im- 
portant modifications. The EH approach contains the minimal number of 
input parameters: current (or pole) quark masses, the string tension a and 
the strong coupling constant a s , and does not contain fitting parameters 
as, e.g., the total subtraction constant in the Hamiltonian. It should be 
useful and attractive to consider expanding of this approach to include 
diquark degrees of freedom with appropriate interactions. The preview of 
this program was done in [14] . It is based on assumption that chiral and 
the short range gluon exchange forces are responsible for the formation 
of ud diquarks in G while the strings are mainly responsible for binding 
constituents in G. In this paper we review and extend application of the 
EH approach to the Jaffe-Wilczek model of pentaquarks. 

In this model inside 0(1540) and other q A q baryons the four quarks are 
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bound into two scalar, singlet isospin diquarks. Diquarks must couple to 
3 C to join in a color singlet hadron. In the quark model five quarks are 
connected by seven strings. In the diquark approximation the short legs 
on this figure shrink to points and the five-quark system effectively reduces 
to the three-body one, studied within the EH approach in [T^l HB]. In total 
there are six flavor symmetric diquark pairs [ud] 2 , [ud][ds] + , [<is] 2 , [ds][si/| + , 
[su] 2 , and [su][i£<i] + combining with the remaining antiquark which give 18 
pentaquark states in 8^ plus 10^. All these states are degenerate in the 
SU(3) F limit. 



2. The EH approach and the results 



The EH for the three constituents has the form 

3/2 



^Eii + fJ+^o + K (1) 

where Hq is the kinetic energy operator, V is the sum of the perturbative 
one-gluon exchange potentials and the string potential String- 

The dynamical masses fa (analogues of the constituent ones) are ex- 
pressed in terms of the current quark masses m a from the condition of the 
minimum of the hadron mass as function of fa 2 -*: 



dM^jm^fa 
dfii 



°. <=E(^+f)+w. ( 2 ) 



Eo(ni) being eigenvalue of the operator Hq+V. Quarks acquire constituent 
masses fa ~ y/a due to the string interaction in (CD). As of today the EH 
in the form of (JH) does not include chiral symmetry breaking effects. A 

2 ) Technically, this is done using the auxiliary field approach to get rid of the square root term in the 
Lagrangian [171118) . Applied to the QCD Lagrangian, this technique yields the EH for hadrons (mesons, 
baryons, pentaquarks) depending on auxiliary fields fa. In practice, these fields are finally treated as 
c- numbers determined from J5J. 
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possible interplay with these effects should be carefully clarified in the 
future. 

The physical mass Mh of a hadron is 

M H = + Y,Ci. (3) 

i 

The (negative) constants C{ have the meaning of the constituent self en- 
ergies and are explicitly expressed in terms of string tension a [19\: 

Q = -—th, (4) 

where 

rj q = 1, r) s = 0.88, rj c = 0.234, rj b = 0.052. (5) 

In Eq. (J3J) 7/ s , rjc, and 775 are the correction factors due to nonvanishing 
current masses of the strange, charm and bottom quarks, respectively. 
The self-energy corrections are due to constituent spin interaction with 
the vacuum background fields and equal zero for any scalar constituent. 

Accuracy of the EH method for the three-quark systems is ~ 100 MeV 
or better JTHl CHj- One can expect the same accuracy for the diquark- 
diquark-(anti)quark system. 

Consider a pentaquark consisting of two identical diquarks with current 
mass m\ u j\ and antiquark with current mass mq (q = d, s, c). In the hyper- 
spherical formalism the wave function ip(p : \) expressed in terms of the 
Jacobi coordinates p and A and can be written in a symbolical shorthand 

as 

^(p,A) = ^^(^)V]( Q )' ( 6 ) 

K 

where Y\jc\ are eigen functions (the hyperspherical harmonics) of 
the angular momentum operator K(Q) on the 6-dimensional sphere: 
K 2 (Q)Y{ K } = —K(K + 4)Yjifi, with K being the grand orbital momen- 
tum. For identical diquarks, like [ud] 2 , the lightest state must have a wave 
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function antisymmetric under diquark space exchange. There are two pos- 
sible pentaquark wave functions antisymmetric under diquark exchange, 
the first one (with lower energy) corresponding to the total orbital mo- 
mentum L = 1, and the second one (with higher energy) corresponding to 
L = 0. For a state with L = 1, l p = 1, l\ = the wave function in the 
lowest hyperspherical approximation K = 1 reads: 



4> = iJ- 5 / 2 xi(fl)tii(«), ui(fl) = J^sinB- Y lm (p), 



(7) 



where i? 2 = p 2 + A 2 . Here one unit of orbital momentum between the 
diquarks is with respect to the p variable whereas the A variable is in 
an (S-state. The Schrodinger equation for Xi{R) written in terms of the 
variable x = ^JjiR, where fi is an arbitrary scale of mass dimension which 
drops off in the final expressions, reads: 



d 2 Xi(x) 
dx 2 



+ 2 



E H bix - — 

x Hx z 



XiW = 0, (8) 



with the boundary condition Xk{%) ~ 0(x 7 ^ 2 ) as x — > and the asymp- 
totic behavior Xi{ x ) ~ Ai((26i) 1 / 3 x) as x — > oo. In Eq. 



(9) 



ai = fl-^/J • y V c (ri, r 2 , r 3 ) • w 2 • dtl, 
h = ■ J String r 2 , r 3 ) • • dfl, 



where 



and 



Vc(^i,r 2 ,r 3 ) = --a s • , (10) 



Vrtring(ri, r 2 , r 3 ) = cr ■ / min (11) 

is proportional to the total length of the strings, i.e., to the sum of the 
distances of (anti)quark or diquarks from the string junction point. In the 
Y-shape, the strings meet at 120° in order to insure the minimum energy. 
This shape moves continuously to a two-legs configuration where the legs 
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meet at an angle larger than 120°. Explicit expression of KtringC^ij ^2? ^3) 
in terms of Jacobi variables is given in [20J. 

The mass of the G + obviously depends on m\ u d] and m s . The current 
masses of the light quarks are relatively well-known: m u ^ ~ 0, m s ~ 170 
MeV. The only other parameter of strong interactions is the effective mass 
of the diquark my. In principle, this mass could be computed dynami- 
cally. Instead, one can tune m\ u j\ (as well as m\ us ] and mu s ]) to obtain the 
baryon masses in the quark-diquark approximation. We shall comment 
on this point later on. 

In what follows, we use a = 0.15 GeV , and explicitly include the 
Coulomb-like interaction between quark and diquarks with a s = 0.39. 

For the pedagogy, let us first assume myi = 0. This assumption leads 
to the lowest uuddd and uudds pentaquarks. If the current diquark masses 
vanish, then the [ud] 2 d pentaquark is dynamically exactly analogous to 
the J P = \ nucleon resonance and [ud] 2 s pentaquark is an analogue of 
the J P — \ A hyperon, with one important exception. The masses of 
P-wave baryons calculated using the EH method acquire the (negative) 
contribution 3C q for J P = \ nucleons and 2C q + C s for the J p = \ 
hyperons. These contributions are due to the interaction of constituent 
spins with the vacuum chromomagnetic field. Using the results of Table 1 
below we get the mass of the P-wave nucleon resonance with the orbital 
p excitation 1600 MeV and the mass of A hyperon 1600 MeV that within 
100 MeV agrees with the known P-wave N and A resonances. 

However, the above discussion shows that the self-energies C\ U( j\ equal 
zero for the scalar diquarks. This means that introducing any scalar con- 
stituent increases the pentaquark energy (relative to the N and A P-wave 
resonances) by 2\C q \ ~ 200 — 300 MeV. Therefore prior any calculation we 
can put the lower bound for the pentaquark in the Jaffe-Wilczek approx- 
imation, M(0) > 2 GeV. 
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The numerical calculation for mu^i = yields the mass of [ud] 2 s pen- 
taquark ~ 2100 MeV (see Table HI). The similar calculations yield the 
mass of [ud] 2 c pentaquark ~ 3250 MeV (for m c = 1.4 GeV) and [ud] 2 b 
pentaquark ~ 6509 MeV (for rrib = 4.8 GeV) [21]. For illustration of ac- 
curacy of the auxiliary field (AF) formalism in Table [T] are also shown 
the masses of [ud] 2 d and [ud] 2 d pentaquarks calculated using the spinless 
Salpeter equation (SSE): 



i=i 

M = M > 



" i=i ( VpI + m l 



where V is the same as in Eq. (CD), Mq is the eigenvalue of Hs, 

( yjpfudi ^ m \ud] } ' (y^P|+~^l} are ^ ne avera g e kinetic energies of di- 
quarks and an antiquark and j]i are the correction factors given in (Ed- 
The numerical algorithm to solve the three-body SSE is based on an ex- 
pansion of the wave function in terms of harmonic oscillator functions with 
different sizes [22] • In fact to apply this techniques to the three-body SSE 
we need to use an approximation of the three-body potential Turing by a 
sum of the two- and one-body potentials, see |23j. This approximation, 
however, introduces the marginal correction to the energy eigenvalues. The 
quantities and fi q denote either the constituent masses calculated in 

the AF formalism using Eq. (JH) or ^ \Jp\ u d] + m \ud] } > \\[p\^~ m f) f° un d 
from the solution of SSE. It is seen from Table 1 that these quantities agree 
with accuracy better than 5%. The pentaquark masses calculated by the 
two methods differ by 100 MeV for ([ud] 2 s) and 160 MeV for [udfd. The 
approximation of String mentioned above introduces the correction to the 
energy eigenvalues < 30 MeV, so we conclude that the results obtained 
using the AF formalism and the SSE agree within ~ 5%, i.e., the accuracy 
of the AF results for pentaquarks is the same as for the qq system (see, 
e.g., |21). 
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If we withdraw an assumption m\ u d\ = 0, then a possible way to estimate 
the current diquark masses is to tune miy, mr ns i and mu s ] from the fit to 
the nucleon and hyperon masses (in the quark-diquark approximation). 
In this way one naturally obtains larger pentaquark masses. We have 
performed such the calculations using the SSE. We briefly investigated 
the sensitivity of the pentaquark mass predictions to the choice of a, the 
strange quark mass m s and diquark masses miy and found M([ud] 2 d) in 
the range 2.2-2.4 GeV, M{[ud] 2 s) ~ 2.4 GeV and M{[us] 2 d) ~ 2.5 GeV. 

Increasing a s up to 0.6 (the value used in the Capstick-Isgur model 
|25j) decreases the [ud] 2 s mass by ~ 120 MeV (see Table ED- We have 
briefly investigated the effect of the hyperfine interaction due to the a 
meson exchange between diquarks and strange antiquark and found that 
it lowers the G + energy by ~ 180 MeV for g^/An ~ 1. As the result we 
obtain the lower bound of [ud] 2 s pentaquark M([ud] 2 s) = 1740 MeV (for 
f^[ud\ = 0) which is still ~ 200 MeV above the experimental value. 

3. Conclusions 

We therefore conclude that the string dynamics alone in its simplified 
form predicts too high masses of pentaquarks. This may indicate on a 
large role of the chiral symmetry breaking effects in light pentaquark sys- 
tems. An "extremal" approach of chiral soliton model totally neglects 
the confinement effects and concentrates on the pure chiral properties of 
baryons. Therefore the existence of O, if confirmed, provides an unique 
possibility to clarify the interplay between the quark and chiral degrees of 
freedom in light baryons. 

This work was supported by RFBR grants No. 03-02-17345, 04-02- 
17263, the grant for leading scientific schools No. 1774.2003.2. The NATO 
is also greatly acknowledged for the grant No. PST.CLG. 978710. 
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Tables 



Table 1. The pentaquark masses in the quark-diquark-diquark approxi- 
mations. The masses of [ud] 2 q states (q = d, s) for J p = ^ + pentaquarks 
are shown in units of GeV. 











M 


[ud] 2 s 


1 + 
2 


AF 


0.482 


0.458 


2.171 






SSE 


0.463 


0.468 


2.070 


[ud] 2 d 


1 + 
2 


AF 


0.476 


0.415 


2.091 






SSE 


0.469 


0.379 


1.934 



Table 2. The [ud] 2 s mass calculated using the SSE as a function of a s 
with and without the Goldstone boson exchange (GBE) in units of GeV. 





0.39 


0.50 


0.60 


M (with GBE) 


1.893 


1.814 


1.737 


M (without GBE) 


2.069 


2.007 


1.949 
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